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We consider a class of warped higher dimensional brane models with topology A4 
where E is a D2 dimensional manifold. Two branes of codimension one are embedded in such a bulk 
' space-time and sit at the orbifold fixed points. We concentrate on the case where an exponential 

I warp factor (depending on the distance along the orbifold) accompanies the Minkowski M and the 

04 , internal space E line elements. We evaluate the moduli effective potential induced by bulk scalar 

fields in these models, and we show that generically this can stabilize the size of the extra dimensions. 
As an application, we consider a scenario where supersymmetry is broken not far below the cutoff 
■ scale, and the hierarchy between the electroweak and the effective Planck scales is generated by a 

combination of redshift and large volume effects. The latter is efficient due to the shrinking of E at 
the negative tension brane, where matter is placed. In this case, we find that the effective potential 
can stabilize the size of the extra dimensions (and the hierarchy) without fine tuning, provided that 
^>f-^ I the internal space E is flat. 
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I. INTRODUCTION 



It is a widely accepted idea that our space-time might have more than four dimensions. This very old proposal 
[l^2], revived in [3], is currently being considered in several contexts from particle physics to cosmology. 
' , From a more fundamental perspective, string theory provides a strong motivation for considering higher dimen- 
sional scenarios. In a suggestive model [4,5], Hofava and Witten considered an 11— dimensional supcrgravity theory 
(— I \ compactified on x Cy x 6*^/^2, where M. represents four dimensional Minkowski space and Cy is a Calabi-Yau 
• space, with the size of the Calabi-Yau much smaller than the orbifold one. In such a set-up, the fixed points of the 
. 5^ ^ Z2 symmetry can be viewed as 10— dimensional boundaries which may accommodate a realistic gauge theory. 

: A closely related model has been considered by Randall and Sundrum [6], where two flat 3— branes of opposite 
$^ ' tension are embedded in 5— dimensional Anti-de Sitter space. The extra dimension is compactified on a Z2 orbifold, 
5^ 1 with the branes sitting at the fixed points. This results in a space-time with a non-factorisable geometry which, 
through a redshift effect, induces a four-dimensional effective scale on one of the branes which is much smaller than 
the 5— dimensional cutoff' scale (which is assumed to be Planckian) . This geometric effect could explain the hierarchy 
between the Planck and the electroweak scales. The "radion" mode (which corresponds to the size of the extra 
dimension) plays a crucial role in this scenario, since its vacuum expectation value (VEV) determines the hierarchy. 
In the original RS model, the radion was massless at the classical level. This would cause an unacceptable modification 
of the resulting gravity on the brane [7], and therefore the radion must be stabilized by some mechanism [8]. Moreover, 
for a satisfactory solution to the hierarchy problem, such mechanism should give the radion a suitably large VEV 
without the need of introducing any fine tuning of parameters. 

The problem of radion stabilization is not particular to the RS model, and arises generically in higher dimensional 
Kaluza-Klein theories, including the usual factorisable geometries based on direct products of Minkowski space with 
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homogeneous manifolds. In that context, Candelas and Weinberg reaHzcd that the Casimir energy of matter fields 
or gravity eould induce an effective potential for the radion at the quantum level, offering a natural mechanism for 
stabilizing the size of the extra dimensions [9]. Following this idea, a number of authors have investigated quantum 
stabilization in the context of brane models. In [10-16] the lowest order quantum corrections arising from bulk fields on 
the RS background have been calculated and their relevance to the radius stabilization problem has been investigated 
(sec also [17-22] for discussions of the possible relevance of quantum effects in cosmological brane-world scenarios and 
[23] for finite temperature effects). Initially, such investigations were somewhat discouraging, since it was found [13] 
that in the RS model it is not possible to stabilize the extra dimension and simultaneously solve the hierarchy problem 
by means of the Casimir forces due to bulk gravitons or generic bulk scalars. However, it has recently been shown 
[24] that bulk gauge fields (or any of their supersymmetric relatives) can do the job. These fields induce logarithmic 
contributions to the radion effective potential which are sufficient to stabilize it, generating a large hierarchy of scales 
without fine-tuning. The logarithmic behaviour can be understood, in a 4D holographic description, as the running 
of gauge couplings with the infrared cut-off scale (which corresponds to the electroweak scale) . 

The RS model has opened up a very interesting framework for model building in particle physics, with possible 
cosmological implications. This scenario, however, is just the simplest possibility within a more general class of higher 
dimensional warped geometries which deserve fuller exploration, specially with regard to their possible embedding 
in string theory. In this connection, one expects that more general internal spaces and higher dimensionalities have 
to be considered, and in this paper we shall take a step in this direction. Our aim is twofold. On one hand, the 
consideration of more general spacetimes may provide interesting extensions of the RS mechanism for the geometric 
origin of the hierarchy. On the other, quantum effects in such scenarios can be qualitatively different, providing new 
ways of stabilizing the radion which do not necessarily rely on the peculiar behaviour of bulk gauge fields. 

Generically, we expect that the behaviour of the effective potential for the "moduli" should be qualitatively different 
once we go beyond the RS scenario. This is indicated (even in five dimensions) by the study of generalized warped 
compactifications [14]. Aside from the non-local Casimir interaction between the branes which we mentioned above, 
local terms which are induced by quantum effects may stabilize the moduli when we consider warped brane worlds 
where the bulk is different from AdS. In the RS model both the branes and the bulk space-time are maximally 
symmetric and thus any possible counter-term amounts to a renormalization of the brane tensions. However, this 
is not true in general. An explicit example was worked out in detail in [14], where a class of 5— dimensional brane 
models with power law warp factors were investigated. In this case, the global symmetry which is responsible for 
the masslessness of the moduli at the classical level is shown to be anomalous. Thus, the effective potential develops 
terms which do not scale appropriately under the global symmetry and which therefore act as stabilizers for the 
moduli. Some of the 5 dimensional models considered in Ref. [14] can be obtained by dimensional reduction of 
5 -|- £'2-diniensional models, and in this paper we shall focus on a class of higher dimensional models which includes 
those. 

Specifically, we shall consider spaces with line element given by 

dsfo+i) = e^^'-y^rj^^dx^dx" + E? e^P^^hijdX' dX^ + dy^. (1) 

The coordinates ,r'^ parametrize four dimensional Minkowski space ^A and the coordinates X' cover a _D2 ^dimensional 
compact internal manifold S. We locate two D = {4 + £^2)— dimensional branes at the fixed points of the orbifolded 
dimension labelled by y. Such a metric is found as a solution of a (D -|- 1)— dimensional system of gravity plus 
certain 'matter' fields. Depending on the field content, different warpings can arise. For instance, Gregory showed 
in [25] that a six dimensional global string solution exists for negative cosmological constant, with a{y) = —k\y\ 
and p{y) = constant. Gherghetta and Shaposhnikov [26] constructed the metric solution of a six dimensional local 
string-like defect with aiy) = p(y) = —k\y\. Generalizations of this models with more extra dimensions with S = S^^ 
were considered in in [27,28] using bulk scalar fields with a hedgehog configuration. 

The authors of [29] include Yang-Mills (YM) fields with appropriately chosen gauge group, instead. They find a 
series of solutions classified in terms of the Ricci flatness of the manifold S: when the internal space is Ricci fiat (for 
example a D2— dimensional torus or a Calabi-Yau space), one obtains warp factors which can generically be expressed 
as sums of exponentials if there is no YM flux: in particular, when the bulk cosmological constant is negative, a 
specially simple solution with both warp factors equal to the RS one exists (including the case of higher dimensional 
AdS space). Turning on some YM flux can relax the condition of the Ricci flatness of the internal manifold. In this 
case, they find a solution where along the Minkowski direction the warp is d la RS, whereas along the curved manifold 
it is constant. This gives the interesting combination of a higher dimensional theory which is a hybrid between an 
ordinary Kaluza-Klein theory and the RS model. In passing, we note that the phenomenology of such a scenario has 
been recently considered in [30,31]. 

In the present paper we start considering quantum effects from bulk fields quantized on such space-times beginning 
with the case when the two warp factors are of the RS type. The computation is similar in some respect to the RS 
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one with some additional technical complications. Moreover we allow a more general setting than the one presented 
in [29] in that we consider solutions with internal manifolds E which arc not Ricci flat. In section II we describe how 
to construct these types of models (for the case of Ricci flat internal manifolds we refer the reader to reference [29]) 
and subsequently discuss the moduli approximation in this type of space-times. 

The Kaluza-Klein reduction of the bulk field is considered in section III, where it is shown that, when branes have a 
non trivial topology, the physical mass depends also on the eigenvalues of the Klein-Gordon operator on the internal 
manifold. The one-loop effective action is evaluated in section IV and the result is simply expressed in terms of 
heat-kernel coefficients. This is done using the Mittag-LefHer expansion of the generalized ( function. The result is 
regulated by using dimensional rcgularization and appropriately renormalized in section V. 

In section VI, we propose a scenario where supersymmetry is broken just below the cutoff scale and the hierarchy 
is generated by a combination of redshift and large volume effects. The possibility of stabilizing the moduli (and the 
hierarchy) with the quantum effective potential is discussed in section VII for the different cases of interest and it is 
shown that such a stabilization can in fact take place without fine tuning. The conclusions arc left to the last section. 

The paper is equipped with an appendix in which we obtain the Mittag-Lcfflcr expansion for a positive semidefinite 
operator (This last result does not claim any originality and can possibly be found elsewhere). Appendix B contains 
the details of the dimensional reduction of the gravitational part of the action. 



II. SETUP 



We are interested in the quantum effective action arising from a quantized bulk scalar field on the background 
space-time specified by the metric (1). We consider two branes of codimension one with the topology of Ad'^ x S, 
where the manifold E is taken to be Einstein and compact. The branes sit at the orbifold fixed points and the Z2 
symmetry is imposed on the solutions. 

As we have already mentioned in the introduction, such solutions have been obtained in [29] for Ricci flat internal 
manifolds. However, as also shown in [27,28], more general solutions with E = S^^ can be foimd by introducing 
additional matter content coming from a scalar field with a hedgehog configuration. In the next Section, we show 
how to obtain such a solution. 

A. Model 

The specific space-time that we consider in this paper corresponds to the case when the two warp factors are equal 
and exponential, <j{y) — p{y) — ~k\y\. The model consists of an G invariant non-linear sigma model parametrized by 
a set of bulk scalar fields together with a standard bulk gravity sector, and two boundary-branes. This is described 
by the action 

S = J d'^+'x^/g^l^M'^-'n^n+i) - A - ^m.^^^^^^ - A(<^«t^a _ ^2^| 

- j d^x^g(^D)+ '^+~ j d^^V9iD)- T-. (2) 

Our notation the following. The higher dimensional bulk indices are M, N, . . . and run over /x, i, y; the (4 -I- D2) = D 

dimensional brane indices are A,B,... and run over /x, i; Q^mn^^ ^he bulk metric and Q^ab"^ ^'^^ ^he induced metrics 
on the branes. Finally, r± are the brane tensions, and M is the higher dimensional fundamental Planck mass. 

Let us look more closely at the structure of the scalar fields. The equation of motion for the scalars can be written 
as usual: 

□0" = -A(^". (3) 
The role of non-dynamical auxiliary field A is to impose the constraint 

</,°t^« - = 0. 

Differentiating this constraint twice, we can rewrite Eq. (3) as follows: 
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The previous equation allows hedgehog solutions for for suitable choices of the group G. Moreover, they have a 
constant profile along the orbifold and satisfy 

7" 2 2 

A^^!.« = -LV", and dM(l>"'^ d'^ (j)" = e-'^P (5) 

where i is a 'winding number', and is the laplacian obtained from 7^ . 

The Einstein equations for such hedgehog configurations have been studied in [26-28] , where solutions of the type 
(1) with a{y) — p{y) — —ky have been found, with 

k = Y^-4Mi-DA/(D - !)£», (6) 

where A < 0. In order to obtain the space-time described previously we take two copies of a slice of this D + 1 
dimensional space comprised between and corresponding to the brane locations. The two copies are glued 
together there. Along with the the identification y — y± — > 2y± — y, this gives the topology of an S^/Z2 orbifold in 
the y direction. 

In order for this to be a solution of our model (2), the brane tensions have to satisfy 



T± = ±A^-{D - VjM^-^K/D = ±4{D - l)M^-^k, (7) 

as a result of the junction conditions at the brancs. Besides (6), the Einstein equations in the bulk relate the hedgehog 
parameters to the curvature of the internal manifold E as 

= ^M^-^ (8) 

Here, since S is homogeneous, the dimcnsionlcss constant C is defined through Tl^]^ — C^ij, and T^l]^ is the Ricci 
tensor computed out of 7,j . 

Associated to the sigma model scalars a number of Nambu Goldstone modes will be present. However, we shall 
assume that these couple to matter only through gravity, so that their effects are negligible. 

B. Moduli 

One interesting feature of this ansatz is that the parameter R, describing the volume of S, does not appear in the 
equations of motion even in the case of a curved internal space. Moreover, the positions of both branes are free at 

the classical level. They correspond to flat directions in the action and thus are the relevant degrees of freedom at 
low energies. In the moduli approximation, which we shall follow here, they arc promoted to four dimensional scalar 
fields. 

One crucial difference of these solutions with respect to the RS model is that they are not homogeneous along the 
orbifold, even in the case when S is a torus ^. This is due to the compactness of S. In contrast with the RS model, 
the positions of both branes are physically meaningful. 

However, it is clear that a scaling of R is equivalent to a shift in the positions of the branes y±. Therefore, they 
are not independent. Rather, only two moduli are needed. Since we will use several combinations of the moduli 
along the paper, we summarize them briefly now: a± = e~^^± , the physical radii of S at the branes R± = a±R, the 
corresponding dimensionless values r± = a±kR, and a = e~''^^-~'^+^ = a-/a+. 

In addition to the moduli, the massless sector also contains the graviton zero mode. To take it into account, we 
perturb the background solution (1) as follows 

ds^ = dy"^ + e^-^^^) [g^^{x)dx^'dx'' + R\jdX'dX^] . (9) 

Substituting this metric back into the action (2) we obtain the kinetic term for g coming from the bulk part (see [14]). 
The kinetic terms for the moduli y± come from the boundary terms. A computation analogous to that in [14] gives 



5(4) = -m% J d'x^g I [^l - <^2_] n - 4^-1 [[d^^f - I , 



(10) 



^In this case, the solution corresponds to a toroidal compactification of a higher dimensional AdS space. 
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where = a£ ^ = e ^ ^^^^^ ^j^^ effective four dimensional Planck mass is given by 

^ -v^K'^'M'^-yk , with Wj, = 

's 



-v^R^^M^-^/k , with = / V7c^^'^- 



(11) 



We note that the moduU (p± are Brans-Dicke (BD) fields and in the frame defined by g^,^, the kinetic term for (p+ 
has the wrong sign. Introducing the new variables (p and ijj [14,32], 

= if cosh?/; and ip^ = ip ainhijj, 

the Einstein frame is given by g^i, = f^g^u- In this frame the action (10) takes the form 



'(4) 



Dn 



D2 + 2' 



(12) 



and now the kinetic terms are both positive definite. Moreover, we note that the modulus (p decouples in the limit 
D2 — > 0, as expected, since this case corresponds to the usual RS model, where only one modulus is present. 

We are assuming that the {D dimensional) matter fields x±'^ are localized on each brane and so they couple 
universally to the corresponding induced metrics g'^^ ^ (recall A,B,--- = n,i) 



Here, we have kept the E— zero modes only, and integrated out the X dependence, the S volume factor has been 
absorbed by the four dimensional matter fields x ^nd couplings, and the four dimensional induced metrics are the 
(/X, u) components of the D dimensional ones g"^^ = g^S^"^ . 
A repeated use of the chain rule leads to the interaction of the moduli with matter given by 



^mod-matt ^ / ^4^ 




Do 



£>2 + 2 V 



2C± 



6hnp ■ 



D2 + 2 



'2+2)/2 



T± + D2C± 



(14) 



where and C± are defined according to ^—g± T± = \/—gT±, and \/—g±a^^ C± = \p-g t±- The coupling of 
the moduli to the Lagrangian is entirely due to the dimensions along E being warped, and is a generic prediction of 
models with a nontrivial warp factor for the extra dimensions. In fact, the 'radion' modulus \E' is coupled to matter 
through (T + DiC)±^ which coincides with the trace of the D dimensional energy momentum tensor. Moreover, this 
shows that the the modulus <~p decouples from matter in the RS limit — > 0, as it should. 
Defining the canonical fields 



■ mp (51n<^, 



£'2 + 2 

we obtain the equations of motion for the moduli 



and 



^ ^2+3 



2 V -D2 + 2 Top 



T+ - 2£+ + T_ - 2£_ 



V^2(D2 + 3)(i:'2 + 2)TOp 



,(D2+2)/2 



(15) 



As we explain in Sec. VI, we are interested in the case of a <C 1 in order to have a substantial redshift effect arising 
from the warp factors. Unless otherwise stated, we shall set (a+) = 1, so that, with a good accuracy, a_ ~ a, 
(/? ~ ifj^ ~ 1 and -0 ~ <C 1. 

Thus, from (15) we can read off the couplings to the two types of matter:^ $ couples to the matter at either brane 
X±, with a strength ^ 1/mp. As for the coupling to X- is quite large, of order a~^^^+^^''^/TOp, and to x+ is even 
smaller than Planckian, <^ a^^^+^^/^/TOp. 



In the rest of this article, we will consider only matter located on the negative tension brane. Here we just consider other 
possible forms of matter aX, y = y+ for the sake of generality. 
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III. KALUZA-KLEIN REDUCTION 



Before the evaluation of the one-loop effective action, which we do in Section IV, we turn now to the reduction in 
KK modes of a bulk scalar field living in the space-time described in the previous section. 

The idea is very simple: by performing a Kaluza-Klein reduction of the higher dimensional scalar field theory from 
D + 1 (with D = Di+ D2) down to Di dimensions, we obtain an equivalent lower dimensional theory consisting of an 
infinite number of massive Kaluza-Klein modes. Specifically, the Kaluza-Klcin reduction is performed by expanding 
the higher dimensional scalar field in terms of a complete and orthogonal set of modes and then integrating out the 
dependence on the extra dimensions. The masses turn out to be quantized according to some eigenvalue problem and 
depend on the details of the space-time, the nature of the internal manifold and on the bulk (higher dimensional) 
scalar field. The one-loop effective action can then be evaluated by re-summing the contribution of each one of the 
modes. 

Typically in Kaluza-Klein theory the mass eigenvalues are found explicitly and the subsequent evaluation of the 
sum over the modes does not present particular difficulties. However, in the case of warped space-times the main 
difference is that the orbifold nature of the extra dimension complicates the mass eigenvalues, which are expressed in 
terms of a transcendental equation and thus cannot be found explicitly. 

In the present section we will carry out the first step of the computation, namely the Kaluza-Klein reduction of the 
bulk scalar field. We will consider the most general case of a massive non- minimally coupled scalar field and assume 
that S, a compact manifold. 

The bulk scalar field T{X,x,y) obeys the following equation of motion: 

[-□(D+i)+"^'+?7^(D+i)]Tr = , (16) 

where 7^(0+1) is the higher dimensional curvature and □(£)_!_ is the D'Alcmbertian, both computed from the metric 
(!)• 

Using the explicit expression for the metric tensor, we can disentangle, in equation (16), the dependence on the 
internal manifold from the Minkowskian one. A straightforward calculation gives: 



E? 

+^2 ^ ^g-2p^^(7) _ ^p^y^ ] T = , (17) 
where A(^) is the laplacian related to 7ij, □ is the Di dimensional flat D'Alcmbertian, and 

F{y) = 2r"(y) + t' {yf + D^ivf + D^p'iyf , 

T{y) = D^a{y) + D2p{y) . (18) 

We now expand the field T(.t, X, y) in terms of a complete set of modes carrying a momentum along the orbifold and 
S directions labelled by indexes n and I respectively, 

T(x, X,y) = Y^ vi/j(X)$,,„(a;)Z,,„(y). (19) 

Here, the modes '^i{X) are a complete set of solutions of the Klein- Gordon equation on the manifold S: 

Pe*;(X) ^^[- A(^) + ^i{X) = Xf^iiX) , (20) 

with eigenvalues and degeneracy"^ gi. If we now require $;,„(x) to satisfy the Klein-Gordon equation on the 
Minkowskian factor of the space- time M. with masses mf^, 

[-a + ml„]^i,n{x)=0 , (21) 



^Although wc assume to be citlier positive semidcfiriite or positive definite, tiie label 1 = always refers to the zero 
eigenvalue, i.e., Ao = 0, the existence of this eigenvalue being set by go being or 1. 
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we are left with a radial equation for the modes Zi^niu) of the form 

e^" [-e-^dye^dy +m^- ^F{y) + Xfe-^"] Z^^^ = ml„Zi,n • (22) 

This equation is valid for any warp factors a and p, and can be viewed as an eigenvalue problem for the orbifold 
modes Zi,n and the physical masses m; „. Both of them depend in general on the 'internal' index I. In this paper we 
consider the case of two equal warp factors, with 

P{y) = o-(y) = -k\y\ . 

Defining D = Di + D2 , we can specialize Eq. (22) to this case as 

-e(2-^)''a,e^-a, + m^e^- - ^F{y)e^-] - {ml„ - Xf)Zu,. (23) 

We note that the operator in the l.h.s. does not depend on the internal index I. Accordingly, in this case neither the 
modes nor the combination qf^ = — depend on Z. In other words, the dependence on I and n of the masses 
is factorized for this geometry, 

mln = Ql + >^l (24) 

Therefore, from now on we shall drop this index in Z. On the other hand, Eq. (23) is similar to the one which arises 
in the RS model, and the most general solution can still be written in terms of Bessel functions: 



Z^{y) = e,{y) [a^J. [^e-^) + B^Y. (|e-)] (25) 



k J ■ ^ -\k 

where for notational convenience we have defined 



and 

r^' = '^-D{D + l)^+^. (27) 

The index /3 has been introduced in order to discriminate the two possible cases of T being untwisted (Z„(— y) — Zn{y)) 
or twisted (Z„(— y) = —Zn{y)). Imposing the appropriate boundary conditions, which can be obtained by integrating 
equation (23) across the orbifold fixed points, we find that the eigenvalues g(„ are determined by the transcendental 
equation: 



F«'(a=0. (28) 



^ka 
where 

pf:)t^\ = fY„{az)J^{z) - J„{az)Y„{z) p = twisted , ^29) 
\y:,{az)j^{z) - j„{az)y„{z) /3 = untwisted . ^ ' 

As in the RS model, the combinations of Bessel functions relevant to the untwisted case are given by 

j,{z) = ^Dil - AOM^) + zJliz) , 
y4z) = ^D{l-40Y4z) + zY:,{z) , 

This completes the Kaluza-Klein reduction of the bulk scalar field. 

In the following we will report only on the case of untwisted fields, although the case of twisted fields can be 
obtained at ease with simple modifications of our calculation. 
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IV. QUANTUM EFFECTIVE ACTION 



The one loop effective action F can be expressed as the sum over the contributions of each mode, r;^„ 

The previous expression can be evaluated in a variety of ways (see for instance [9,33]). Dimensional rcgularization of 
the 4— dimensional Minkowski directions to 4 — 2e leads to the following expression for the vacuum energy contribution 
to the effective action 



with 

y-g(s) = -\{A^y^,■'^T{s)Y^g^{ql + A?)- , (31) 



where the prime in the sum assumes that the zero mass mode is excluded (since it does not contribute) and s = — 2 + e. 
The renormalization scale ^ is introduced for dimensional reasons. It is convenient to separate V^'^^ into three 
contributions^ 



where 



y-s(s) = V^{s) + Vrs{s) + K(s) , (32) 



2e °° 



VRsis) = -go^7^{ka/f,r''Tis)^x-'^ , (34) 

^ ' n=l 
4 oo 

with Xn = Qn/ka and yi = Xi/ka. Thus, Vs retains the contributions from the orbifold zero mode (present only 
for an untwisted field), Vrs is the contribution from the S zero mode (which coincides with the potential in the RS 
model), and 14 includes the contribution from mixed states. It is clear from Eqs. (28) and (29) that x„ depends on a 
only. Since A; scales like 1/R, we can factor out the dependence on this modulus, defining dimensionless eigenvalues 
A; = RXi, that do not depend on R. If S is a one-parameter space, then A; cannot depend on any other shape moduli. 
However, hero wc arc interested in the dependence on the breathing mode R only. So, in general, yi = Xi/{kaR) 
depends on the moduli described in Sec. (II B) through 

The first term in (32) Ve results from the KK excitations along the internal manifold. It can be expressed in terms 
of the generalized C function associated to the laplacian Ps defined on S (see Eq. (20)), ^ 

oo 

as)^as\PE) = j29iK''- (36) 
1=1 

using the previous rescaling we can recast Ve as 

Ve{-2 + e) = if^Rf' r(-2 + e)C(-2 + e) , (37) 



*Here we define Ao = 0, so that the existence of such a zero eigenvahie or not is controlled by go. If go 7^ 0, the RS contribution 
comes about explicitly and introduces a divergence which needs to be cancelled by a corresponding contribution coming from 
V*. This cancellation provides a non-trivial check of our evaluation. The case of a strictly positive definite operator, can be 

obtained by putting go to zero. 

^The fact that this term does not depend on the eigenvalues q„ is a consequence of how we have performed the Kaluza-Klein 
reduction. 
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where we redefined the renormaUzation constant /i. The previous expression can be elegantly dealt with by using the 
Mittag-Leffler representation for the C, function, which proves to be a very useful tool to handle the pole structure of 
the ( function, since the residues at the poles are determined by geometrical quantities of S (See for example [34]). 
As shown in appendix A, 



1 



r{s) 



E 

^P=0 



Cry 



S-D2/2+P 



(38) 



where Cp = Cp — go ^p^D-^ji and Cp are the (integrated) Seeley-DeWitt coefficients of the operator Ps on E, p runs 
over the positive half integers and /(s) is an entire function. In fact, the sum (38) runs over half integers, but, since S 
has no internal boundaries, the coefficients Cj/2 are zero. Relation (38) can now be used to regulate Vs, and a simple 
calculation gives 



y5:(s = -2 + e) = - 



327r2i?4 



^-2 + (^^2/2+2 



(39) 



where r2_2 is the constant term in the power series of r(s)C(s) around s = —2 (see Eq. (A17)). 

The term proportional to the RS contribution has been computed in [10,13,16]. Without going into details, we 
write such term as follows: 



where we have introduced 



V(a) 



dzz-^ In 1 



K{az) i„{z) 



(40) 



(41) 



and the coefficients ci and C2 do not depend on a. Here, the coefficient 64 depends on the mass and non-mininal 
coupling of T, and is defined through Eqs. (51,47,45). 

Before entering into the discussion of the higher dimensional contribution due to the mixed KK states V* , we can 
foresee now some of the details of the computation. As mentioned above, the case when E is a torus corresponds to 
a toroidal compactification of a slice of higher dimensional AdS space. Since it is a maximally symmetric space, all 
the geometric invariants are constant, and so proportional to the brane tensions. Thus, the only possible divergence 
that can appear is of the form J d'^x{l + a^). However, regardless of the dimension of E, the contribution from Vrs 
contains a divergence of the form J d'^x(l + a^). Of course, what happens is that aside from the higher dimensional 
divergence, V^, also contains another divergence that cancels the RS one. This feature occurs not only when E is a 
torus. Rather, it is completely general. As we show next and in appendix A, the divergence of K + Vrs proportional 
to / d'^x{l + a^) is always controlled by a geometric invariant related to E (which trivially vanishes for a torus). 
This ensures that if this divergence persists, it is because one can build some operator that behaves like it in this 
background. 

Let us now turn to the evaluation of K . First of all, let us concentrate on the sum 



r{s) J2 91 {xl + yfY 



(42) 



..i=l 



This is not straightforward to compute, however the method developed in [35-37] allows us to perform such a calcu- 
lation. Since, in our case, the evaluation does not present any particular difficulty, we will be brief and address the 
reader to the original references for an introduction to the details of the method. 

The residue theorem permits us to express the sum (42) as a contour integral and an appropriate choice of the 
contour of integration leaves us with 



(43) 



which, by changing variable and by using some known properties of the Bessel functions can be recast as 



E 9i vi^' I (^' - 1)-^ ^ In iMmz)] dz, 



(44) 
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where 



and 



Pi,{z) = Fv{iz) = - [k„{z)i„{az) - K{az)iv{z)] , 



i^z) = zlUz) + -D{1 - 40Iu{z) 
K{z) = zKUz) + ^D{1 - iOKAz). 



(45) 



The integral (43) is considered in [38] for s = 0. In that case, the result can be expressed in closed form. However, 
for s = —2 this does not seem to be possible. 

We can regulate relation (44) using the asymptotic expansions for the Bessel functions. The large z behaviour of 

and can be written as (see e.g. [10]) 



V(^) = 
ku{z) = 



-'e^''\z) 



(46) 



where e('=)(^) = O'^'^-z) is a power series in 1/z beginning with 1. Thus, we can recast the integrand of Eq. (44) in 
the form 



P^iviz) = -^yize^'-^y''Q(^{yiz)e'^'\ayiz) 

TT 



1 



k„{yiz)i„{ayiz) 
iv{yiz)ki,{ayiz) 



(47) 



Up to a constant term, InP^ can be split as 

In [p^imz)] = n\'\z) + nf\z)+nf'\z) 



with 



n\'\z)^\nz+{l~a)yiz , 
nf\z) = \n\Q^^{yiz)Q^>'\ayiz) 



'Hf\z)=\n 



k„{yiz)i„{ayiz) 
iv{yiz)k,,{ayiz) 



(48) 



and correspondingly. 



K(s) = Fi'^s) + vy\s) + y;"^(s). 



(2)/ 



with 



2(47r)- 



:(Wm)' 



-2s 



f^^^J^m'']^ {z'-ir'^^ln[n<r\z)]dz (a = l,2,3). (49) 



The evaluation of Vy\s) is analogous to the one for Ve(s) and, once more, the Mittag-LefHer expansion allows to 
express the result in terms of the heat-kernel coefficients of the operator on S. We find 



Vy\s = -2 + e) = - 



1 



327r2i?4 



1 , 1 1 



C5/2+D2/2- 
1 - al 



kaR 



(50) 
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The second term vj'^\s) can be evaluated^ using the expUcit form of 0^'^ and G^*': 



In 



oo 

E 



{_iy 



dj z- 



for ^; > 1, 



(51) 



the cocfScients dj can be obtained by simply Taylor expanding the logarithm. Using (51) and treating the sum over 
the eigenvalues yi as in the case of Vy, (see App. (A)), we can write V* as 



1 d ( 1 1 

^ + = 32^ ^^^^'^ E Y{f/2j I tC'2+i?./2-i/2 - 9oSi,j] - + n_2+,/2j {{kaRy + {-kRy) . (52) 



t r(j/2) 



The third term in vj^^ (s) is finite by construction, and we can put safely s 



1 



647r2i?4 ^ 



(3) 



(z)dz 



(53) 



Combining the previous results, we obtain the unrenormalized Casimir energy: 
1 



yreg 



327r2ij4 



^ {{kR.y^{-kR^y\\ 

.7=-l 



7^- + (/?j- 



<?o(fci?)Nci 



4 

a C2 



50C?4'54,i)^ (m-R)^^ 



(fc/M)- 



2Y,gi\\Vi{a,R.) 



where 



r (l/20F)C5/2+D,/2 for J = -1 
= <^ (3/2) C2+z,,/2 for J = 

[ - («^j/r(j/2)) C2_j/2+r>2/2 otherwise, 

and we understand that the Seeley-DeWitt coefficients Cj are zero if i < 0, 

r (l/2V^)17_5/2 for j = -1 
7j = <^ (3/2)r!_2 for j = 

[ -((ij/r(j/2))^/2-2 otherwise. 



and 



Vi(a,i?_) 



/OO 



l)ln 1- 



kujyiz) ii,{yiaz) 
kviviaz) iu{yiz) 



(54) 



(55) 



(56) 



(57) 



Equation (54) shows that as wc advanced above, the lower dimensional divergence coming from the RS contribution 
Vrs is always cancelled, independently of the structure of the internal manifold S. On the other hand, the contribution 
from the KK modes along S only (lower dimensional, as well) may give a divergence corresponding to j = 0. This is 
controlled by the Seclcy-DeWitt coefficient C]j^/2+2i a-^d gives 1/2 of the resulting 3/2 factor in f3o, the rest coming 
from the mixed states in V^. In particular, if D2 is odd, then there is no such divergence (if S is boundaryless) , in 
accordance with the absence of any operator that behaves as f d'^x in the background, in this case. 

To conclude this section, we briefly comment on the differences appearing when we consider a twisted bulk field. 
First, since there is no orbifold zero mode, its contribution Vs is not present. One can show that the asymptotic 
behaviour of the function P,^ differs in two powers of the argument, originating a change of sign in the contribution 
to /3o and 70 from fJ^'. Of course, the dj coefficients also change, and can be read from [10,13]. In brief, one needs 
to change the dj by the corresponding one, and the 3/2 factor in /3o and 70 by —1/2. Also, we haven't included any 
brane mass terms or kinetic terms, relevant for the untwisted case only (aside from the ones arising from the coupling 
to curvature). In principle, these can be different on each brane. This changes our result in that we would have 
different coefficients, , for the r± series. 



^Strictly speaking wc arc using an asymptotic expansion and tlicrcforc the equality sign is not exact. However, tfie approxi- 
mation we are making is reasonable because tlie integration range vary from 1 to 00 and the argument of Q*-*' and G'*'-' is large 
in the region J? <C 1 and aR <C 1. 
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V. RENORMALIZATION 



In the previous section, we have computed the unrenormaUzed Casimir energy (54) using dimensional regularization. 
This allows us to isolate the divergent terms, of the form 



-^div 



J=-l 



(58) 



with (3j given by (55). A finite number of divergences appear because we have computed the one loop contribution 
to the effective potential. 

It is well known (see e.g. [37]) that the divergences present in the effective action are given by the Seeley-DeWitt 
coefficient C(£)_|_i)/2 related to the operator in (16) on our D + 1 dimensional background space-time. Since this 
has boundaries, nonzero boundary terms are present for any dimension. Moreover, since the extrinsic curvature is 
constant in the space-time we are considering, several powers of the intrinsic curvature of the boundaries are present. 
Finally, it is easily shown that once any possible bulk term is evaluated on the background solution, it can be recast 
as boundary term for this specific solution.^ 

So, we shall consider boundary term of the form 



N 



(59) 



where N = 0,1,2,... and TZ(d)± denotes the (intrinsic) curvature computed from the induced metrics on the 
branes g[D)±- Using the explicit expression for the metric tensor (1), a simple calculation shows the previous term 
generates a contribution proportional to Jl^^^-^N coming from the brane at y+, and a contribution of the form 
j^D2-2N ^A+D2-2e~2N fj-Qj-^ other brane. Then, it is clear that all the divergences in (58) can be dealt with oper- 
ators of the form (59). Specifically, we can take the following expression as the counter-term needed to renormalize 
the effective action: 



1 



327r2e 



(60) 



The index j here runs over the integers comprised between —1 and 4 -|- D2, and are renormahzation constants. 
We recall that, from (58) and (55), the divergences occur for j even only if is even, and for j odd when D2 odd. 

In the process of subtracting the counter-terms, finite contributions to the vacuum energy with a logarithmic 
dependence on the moduli are generated. The renormalized expression can be written as 



V{R±) 



1 



327r2i?4 



^ < (/3j - 50 di 5i,j) \{kR-y \n.{kR-f + {-kR+y \n{kR+f 



L j=-i 



+ (7^ - /3,- In {k/iif) [{kR_y + {-kR+y] I 

+ ga{kRf {ci + a^C2 - 2a^V{a)] + 2j29i>^tVi{a,R-) 



1=1 



(61) 



A few remarks are now in order. First of all, note that we recast the result in order to isolate the fj, dependent 
terms. Such terms are not computable from our effective theory, rather they have to be fixed by imposing a set 
of renormahzation conditions. Secondly, notice that the result is valid for D2 even as well as for D2 odd, and the 
heat-kernel coefficients automatically take this into account. Also, it should be noted that the flat space limit A; ^ 
is well defined. One has to take into account that R+ = i?e~'^^+ and a similar expression for Then in the limit 



''For instance, / d^^+^-'a; ^5(d+i)A = ^_|_ / dPx ^g(B)±<y±, with cr± = ^2A/Dk. 
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k 0, keeping the proper distance between two branes {y- — y+) fixed, all terms which would be singular in k cancel 
each other in the first and second lines of Eq. (61). The rest of the terms are well behaved. The resulting expression 
is proportional to (y_ — y+)/R^, which has a clear interpretation in terms of the Casimir energy. 

An important remark concerns the divergence proportional to i?^ + i?l . This is the divergence present in the RS 
contribution [10,13]. For D2 odd, it is not reproduced by any of the counter-terms in (59). However, this is not a 
problem because such divergence is cancelled by the corresponding one coming from (52) for j = 4. 



VI. PHYSICAL SCALES AND THE EW/PLANCK HIERARCHY 



In this section we propose a scenario where supersymmetry is broken at a scale ?7susy not far below the cutoff 
scale M, and the hierarchy between the electroweak and the effective Planck scales is generated by a combination of 
redshift and large volume effects. Also, we discuss the range of possible values for the dynamical (the moduli R±) 
and the fixed scales (the cutoff M and the SUSY breaking scale jysusy)- 

From Eqs. (12) and (11), we see that the relation between the four dimensional effective Planck mass and the 
higher dimensional one (in the four dimensional effective theory using the Einstein frame metric g^^i/) is 



(62) 



We shall assume that the masses of particles (located at y — y^) are somewhat below the cutoff M. In the four 
dimensional theory, these masses are redshifted down to ~ aM. Then, the EW/Planck hierarchy is given by 



k 



10 



-32 



(63) 



Thus, the EW/Planck hierarchy h is explained in this model due to a combination of redshift [6] and large volume [3] 
effects (even though the branes are of codimension 1). The crucial ingredient in order for the large volume effect to 
be efficient (aside from having a long orbifold), is that the additional extra space E exponentially grows as one moves 
away from the negative tension brane (see Fig. 1). In this way, matter is allowed to propagate along a small S, of 
size _R_, whereas gravity is diluted since it propagates through a much larger E, of effective size Since the gauge 
interactions must not be diluted by an analogous effect, we have to assume that the compactification scale on the 
negative tension brane l/i?_ is close to the fundamental cutoff M . ^ 

Our model solves the hierarchy problem in a fashion very similar to the models considered in [39,40], with two 
concentric branes embedded in a noncompact bulk. In this references, the hierarchy and the positions of the branes 
are naturally stabilized by a generalization of the Goldberger and Wise mechanism [8] (see also [41]). 




FIG. 1. Matter can propagate along the additional extra space E of size R-, but gravity samples a much bigger space. 



Keeping only the E-zero mode in the action for a D dimensional Yang-Mills field Fab aX y = y- with coupling constant 
dliD) ~ one obtains / d°x._Jg^ — FABFcDg'^^^_ g'^^^_ ^ / d'^x^g — F^.^Fp^g'"'' g""^ , where we used 

that g^i, ~ g^i,. Thus, the four dimensional YM coupling is identified as 5^^^^ — g'^^^^/R^^ ~ l/(Mii_)^^. 
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Let us now examine the constraints that we have on the moduh and the physical values they can take. First of all, 
we are thinking of an inter-branc distance d — |y_ — y+j somewhat larger than the inverse curvature scale 1/fc of the 
bulk, in order to have a substantial redshift factor a = e^^'^ . On the other hand, the smallest physical length scale 
is given by the size of S at the negative tension brane, i?_ . This cannot be smaller than the fundamental length of 
the theory though, as argued in the previous paragraph, it should be close to it. There is a tighter technical 
restriction which we shall use so that the result for the potential (61) is valid. This, is organized as a power series 
in r± = kR±, and can be trusted only when 1/R+ a bit larger than the curvature scale and 1/R- a factor a above 
(recall that the ratio R-/R+ coincides with the redshift factor a). Incidentally, we remark that this corresponds to 
the physical situation where the size of the internal manifold E is everywhere smaller than the inter-branc distance 
~ 1/fc. ^ So, we must assume a separation between the fundamental cutoff M and the curvature scale k at least of 
order a. This leads to the following scenario. 

Consider a supersymmetric theory where the SUSY breaking scale is given by r/suSY-Then, the bulk cosmological 
constant A ~ k'^M^~^ is expected to be proportional to jyfJ^V' which leads to 



Even if SUSY is broken not far below the cut-off scale, this may lead to a curvature scale k many orders of magnitude 
below M, due to the large exponent in (64). If the moduli R± are stabilized near the values i?+ ^ 1/fc and R- ^ 1/M, 
then a ~ k/M and from (63), the hierarchy is given by 



Note that the required hierarchy is obtained with r/susY within one order of magnitude of the cut-off M for D — 10, 
and less than 3 orders of magnitude below M for D = 5. 

This shows how the problem of the stabilization of a large hierarchy works in this model. Having introduced a 
small separation between the SUSY breaking and the cutoff scales, we obtain a stable very flat warped space-time, 
k <C M. If the potential (61) can stabilize the moduli R± near the values, i?+ ~ 1/fc and i?_ ^ 1/Af, then the 
effective Planck mass is very large as compared to the EW scale. Whether or not the effective potential (61) can do 
this job is addressed in next section. 

Let us discuss the physical scales in the model in some detail. As illustrated in Fig. 1. the branes are of codimcnsion 
1, so that matter (residing on the negative tension brane, at y = y-) can propagate through a physical extra 
dimensional space of size R-. The mass scales on this brane are redshifted by a factor a, thus the mass of the 
first KK excitations of matter fields is 1/R. Then, from collider physics, we have to set the compactification scale 
l/R>TeV, at least. 

In contrast, gravity propagates through the whole bulk space, and its KK spectrum is analogous to that obtained 
in Sec. Ill for a scalar field. In particular, there are three kinds of modes, excited along the orbifold only, along E 
only or along both, as (24) shows. The masses of the first graviton KK modes along E are of order 1/R. However, 
the modes winding along the orbifold only (the E zero mode) have masses given by morb ^ o,k, as in the RS model 
(the curvature scale times the redshift factor). In the approximation of everywhere small E that we are considering, 
kR± <C 1, this means that these modes are a factor a lighter than the modes propagating along E. 

This allows us to assume the SUSY breaking scale t7susy and the cutoff M are such that fc ^ TeU, obtaining quite 
small masses for the graviton orbifold KK modes Worb ~ aTeV. Below, we show that such a small value does not 
conflict with observations, since the coupling of these modes to matter is very suppressed. This is consistent with 
the assumption made above that in the higher dimensional theory the masses of matter fields are near the cutoff M, 
since they are redshifted to aM ~ TeV, which compatible with the electroweak scale. Also, since we are considering 
the limit of everywhere small internal space kR± < 1, setting fc TeV implies that the masses of matter KK fields 
is large enough, 1/R > TeV. 

Thus, from the point of view of the 4 dimensional effective theory, KK modes from the matter fields appear at 
1/R ^ TeV. Since the curvature scale fc of the bulk is close to 1/R, this coincides with the scale where gravity 
becomes higher dimensional. 

In summary, we are lead to consider distribution of scales illustrated in Fig. 2. We set the cutoff M and the SUSY 
breaking scale tjsvsy ^ M such that the curvature scale of the bulk is fc TeV. We assume that some mechanism 



®This means that in a certain range of energy the model is effectively 5 dimensional. In appendix B, we derive the form of 
the dimensionally reduced theory down to 5 dimensions. 




(64) 




(65) 
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can stabilize i?_ near the fundamental length 1/M and i?+ ~ 1/fc. As a consequence, the masses of the graviton KK 
modes along the orbifold are nioib ~ aTeV, and for the modes along S are ~ TeV. 



E 

M > r/suSY a'^TeV 

1/R- 



kr^TeV 



roorb ~ ak ^ aTeV 



FIG. 2. The mass morb of the first KK excitations along the orbifold is much smaller than the mass of the modes excited 
along the internal manifold. There is essentially the same hierarchy between the fundamental cutoff M and the scales that 
determine the solution, k and R. 



It remains to be seen that, indeed, the graviton KK modes along the orbifold are unobservable, in spite of their 
relatively small masses morb ~ ak. We see from Eqs. (25,26) that the (unnormalized) wave function of the KK modes 
grows exponentially as e^'^'^'/^, signaling that the more warped the extra dimensions are, the more localized on the 
negative tension brane these modes are. This implies [30] that the coupling of these graviton KK modes is amplified 
with respect to that of the zero mode 1/mp) by a factor a~^^^ — h^^^'^'^^~^\ Then, they are much more weakly 
coupled to matter ~ l/{10^TeV) than in the RS model 1/TeV). Thus, in spite of their relatively small mass, 
these KK gravitons cannot be seen individually in accelerators. Moreover, since they are associated with only one 
ofF-the-brane dimension, they do not have as large a multiplicity as in the usual large volume mechanism [3], and so 
they do not significantly cool stars. The total rate of emission of any of such gravitons at a given energy E < TeV 
can be estimated as the coupling squared times the number of states with masses lighter than E [3,42], 

( ^ V ^ ^ h(D-4)/iD-l)^ ^_ 

\a'^l'^mp) morb TeV TeV^' 

which is very small for the energies available inside stars. 



VII. STABILIZATION 



The result we have obtained so far (61) is a potential V{R±) for the two moduli describing the background. Using 
r± = kR± , it can be cast as 



^^"■^^ = ~32iR^ + V.{r.)+v{r+,r.)] , 



(66) 



'^'^ Another interesting possibility consists of setting morb ~ TeV so that k ~ a~^TeV and M ~ a~^TeV. This could be 
realized in a scenario with the SUSY breaking scale rjsusY ~ 1/Ji ~ k and the masses of particles of order k, from the [D + 1) 
viewpoint. In this scenario, the EW/Planck hierarchy is given by ~ a^k^/m%. If the moduli arc stabilized so that R+ ^ 1/fc 
and R- > 1/M, then h ~ (fc/M)^^^, thus needing less separation between r;susY and M in order to explain same hierarchy 
h. Moreover, one can see that the potential (61) generates masses for the moduli larger than in the scenario presented so far. 
However, the bulk cosmological constant A would be much larger than r^^^V- 
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where v{R, r) contains the 'non-local' part, and 

■DC 

V±{r±) = Yl (Tl)' {^J-^i + - 9od4Sj,4) ri Inr^ - afri] , (67) 
j=-i 

Here, the coefRcients are understood to be finite renormalization constants, and are nonzero when the corre- 
sponding logarithmic term is nonzero. This is dictatc;(l by ft-j being zero or not {i.e., whether or not such a term is 
divergent), with the sole exception of j = 4. If /?4 = and the laplacian (see Eq. (20)) has one zero eigenvalue, 
go — 1, the logarithmic terms corresponding to j = 4 are not associated to any divergence of the effective action, and 
= 0. This situation arises, for example, when E is a torus. 

Note that the sum goes from — 1 to oo and we recall that from Eq. (55), all the Pj with j > 4 + D2 vanish identically. 
Thus, the term (ijrL appears with j running from —1 to D2 + 4, and the same holds for the terms with (there 
are a finite number of divergent terms) . 

One interesting feature of the effective potential (66) in both cases with D2 even and odd is that the two leading 
terms in the small r± limit (corresponding to j = — 1 , 0) do not depend on the mass m nor the non- minimal coupling 
constant ^. This means that if we consider equal number of fermionic and bosonic degrees of freedom, these terms 
cancel identically even with non supersymmetric masses. From now on, we will focus on this case, one motivation 
being that the models considered here arise mainly in string theories, and the field content of the effective theories 
indeed contain equal number of bosonic and fermionic degrees of freedom. The only change is that the sum in Eq. 
(67) will begin at j = 1 instead of j = — 1. As mentioned above, the effective potential contains a finite number of 
renormalization parameters af. Their values are not computable from our effective theory. Rather, we shall fix them 
by requiring some renormalization conditions, which determine the values for the moduli as well. Since the moduli 
must be stabilized, we demand 

dr+V{r±) = dr_V{r±) = , (68) 
and in order to match the observed value of the effective four dimensional cosmological constant, we shall impose 

V(r±)|„,i„~ 10-^22^1,. (69) 

We are interested in the limit when the size of S is everywhere smaller than the orbifold size, r+ < 1 and r_ ^ 1. 

One can show that in this limit the non-local term v{r±) is exponentially suppressed, and we can approximate the 
potential by the 'local' terms V±{r±). Moreover, since we consider only the positive powers of r± in V±, the potential 
at the minimum is dominated by r+. Then, conditions (68) and (69) reduce to 

Vl{r+)=VL{r-) = 0, and R-^V+{r+)Un ^ lO-^^^m%. (70) 

To investigate whether this potential can stabilize the moduli, we consider separately the cases with fiat and curved 
S. 



^^For instance, consider a the six dimensionaJ example, with Yl = S^. As described in more detail in Subsect. VII A, the 
generalized zeta function is related to the Riemann zeta function. In this case the we can easily work out the asymptotic 
behaviour of the nonlocal contribution due to the mixed KK states V((o, R-) defined in (57). If we keep the first term in the 
asymptotic expansion of the Bessel functions (46) 

g-2(l-o)j/i^^ 

afyfUs (e-'('-)y') + 6 (1 - a)yi Lu [e-'^'-^^^') + 3 U, (e-^(i-">«) } . 

Taking only the first term in the series of the poly-logarithms for small arguments Li„(«) « z, recalling that yi = \i/r- with 
M = I, gi = 2 and summing over Z = 1, 2, . . . , we find to leading order 

00 

\4^,/ r, \ -1-2 -2(l-a)/r_ 

Thus, this contribution is safely negligible in the limit of small internal space size r_ <C 1. 



Vi(o,i?-)~ j dz{z^ -z) In (1- 
= _1 I U (1- 
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A. Flat E 



This case corresponds to a toroidal compactification of a 4 + + 1 dimensional RS model (with two codimension 

on(^ brancs). In this case, all the divergences have the same form, because all geometric invariants are constant and 
thus proportional to the brane tensions. Thus, there will appear a logarithmic term in the (4 + D2)-th power of r±. 
As can also be derived from Eqns. (61), (55) and (56), setting Cj = for all i ^ , and go = 1, there is another 
logarithmic term corresponding to j = 4. 

Thus, the expression for the potential reduces to 



V±{r±) w I T 7ir± + l2r± T ld.r\. + (74 - di lnr±) r± + . . . 
+ (Tl)^+^^/34+z,.ri+^^ Inr^ - a^^jt^^ +...}. 



(71) 



To illustrate better how the stabilization mechanism works in these cases, we shall discuss in more detail the six 
dimensional example with E = 5^. 

The laplacian (20) on this flat manifold is Ps = dg/R"^, and its generalized zeta function (36, A2) is related to the 
Riemann zeta function through 

as\dl) = 2U2s). 

The pole structure of Cr{2s) is easily found and one immediately identifies 



and 



7j 



-4^5/3 




-2/945 

3C^(-4) 

4C;,(-2)rf2 

MjCrU - 4) 
(i-4)(i-2) 



for j = 5, 
otherwise, 



for j = — 1 . 
for j = 0, 
for j = 2, 
for j = 4, 

otherwise. 



(72) 



(73) 



Prom (71), the potential is of the form 

V±{r±) « { T lir± + 72r± T 73r-± + (74 - di Inr^) T p5r± Inr^ 



...}. 



(74) 



As we mentioned above, the renormalization constants af arise from a finite renormalization 5t± of the brane tensions, 



so that af ~ 2TrST±/k^. The size of St± is expected to be set by the SUSY breaking scale ?7susy so that af are large 
in principle. Then, the main contributions to this potential arise from the fifth and the first powers. The extremum 
condition for the r_ modulus can be well approximated by 



IOtt r_ 



Setting the natural value 6t- ~ JyluSY' obtain 



(— ) 

\VSVSY / 



1/2 



M' 



so that indeed R- is stabilized just above the fundamental scale 1/M without fine tuning. 

As for r+, we have two conditions for just one variable, St+. The idea is to use the renormalization constant St+ 
in order to satisfy V+lmin — 0, and then using this value in = 0, the r+ is determined. In order to be consistent, 
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we should obtain r+ < 1. In such a case, wc can foresee from Eq. (74) that if a'^ has to compensate for the potential 
at the minimum, it has to be of order one. But this means that 6t+ is fine tuned to a value ~ instead of JyluSY- 
Imposing explicitly these conditions, we obtain 



DTT ri 



and 



471 1 
r+ ~ 1. 

372 

We can easily check that for a twisted field this ratio is ~ 0.6, in agreement with the assumption wc made above. 
For the untwisted case, this ratio depends on the boundary and bulk masses, so it can be made small gcnerically. In 
conclusion, besides the fine tuning needed in order to match the four dimensional cosmological constant, no tuning is 
needed for the Planck/EW hierarchy in this case. 

A simple computation gives the mass that this potential induces for the canonical moduli $ and ^ of section II B 

where we used 1/7? ^ k ^ TeV and Eq. (65). The mass of <I> is of the order of the inverse millimeter, which is 
large enough in order not to cause deviations from Newton's law at short distances. Since, as shown in section II B, 
the coupling of $ to matter is suppressed by a Planckian factor, its effects in accelerators are negligible as well. On 
the other hand, the mass for the modulus * is of lOKeV size. From (15), its coupling to matter is suppressed as 
l/{10^TeV). 

Let us briefly discuss the stabilization when wc consider a liiglic;r dimensional flat S. As mentioned above, with 
more flat dimensions, the renormalization constants related to the brane tensions appear with higher powers of 
r±. The only change with respect to the case above is that the condition V!_ =0 now reads 

r_ 

and assuming a natural value for (5t_ given by ?7|^gYi obtain again r_ ~ (A:/M)(M/?7susy)^^^- Thus, for any 
dimension D the modulus i?_ is stabilized without fine tuning near l/M. 

As for the modulus we expect the potential (71) to stabilize it near k once the fine tuning of 6t+ ~ k^ needed 
for the cosmological constant is performed. 

We can compiite the masses for the moduli for an arbitrary number of flat internal dimensions. We find that the 
mass for the $ is always millimetric, whereas ~ aTeV increases with D2, ranging from lOKeV for D2 = 1 to 
lOOMeV for D2 = 6. The coupling of * to matter, of strength (see Eq. (15)) 



1 

is comprised between ~ l/{10'^TeV) for D2 = 1 and ^ l/(100TeV^) for large D2. This guarantees that it hasn't been 
produced at colliders, or has any eS'ect in star cooling. 



B. Curved S 



When E is not flat, besides the divergences proportional to brane tensions terms (giving rise to the power in 
V±), the potential has more divergences. For instance, there can appear divergences proportional to curvature terms, 
which give rise to the powers Accordingly, terms with fewer powers of r± are due to higher powers of the 

curvature, and in general the effective potential takes the form 

00 D 

V±(r±) = ^(Tl)^ {ijri + (Pjri - god^S,,^) Inr^} -^afrl (77) 
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As in the previous case for the brane tensions, the size of the renormaHzation constants in front of these operators 
are expected to be of order the cutoff scale M (or ?7susy)- Finite renormahzation terms of boundary operators behave 
as, 

and we conclude that the dimensionless renormahzation constants in (77) are large, ~ {M/kY » 1. Thus, these 
terms are a series in MR± > 1 rather than in kR± < 1, the dominant terms being with the highest powers, that 
is, the brane tension and the curvature terms. As a first approximation, we can neglect the remaining terms, and 
minimum condition for R- is reached naturally for _R_ ^ ^/M, which is what we need (see Fig. 2). 

However, we see that in order to obtain ~ 1/k, we need to tune the ratio oi ao and aD-2- Besides, the tuning 
corresponding to the cosmological constant is still needed. 

In principle, we could consider the case when the heat kernel coefficient Ci(Pe) is zero, which can happen for 
some value of the non-minimal coupling ^. We see from (55) that in this case there is no divergence in the potential 
corresponding to the curvature terms. Then, assuming that the next nonzero coefficient is C2, the two powers that 
dominate the potential are (MR±)^ and (MR±)^~'^ . However, in order to stabilize i?+ near 1/k, again we have to 
do one fine tuning. We can say that in general, the presence of any other divergence, besides the brane tension, spoils 
the efficiency of the potential in stabilizing the moduli at well separated scales. 

We conclude that, for curved S the potential can naturally stabilize the moduli but without a large hierarchy. 

VIII. CONCLUSIONS 

In this article we have investigated the role of quantum effects arising from bulk scalar fields in higher dimensional 
brane models. Specifically, we have considered a class of warped brane models whose topology is x S x S^/Z2, 
where E is a D2 dimensional one-parameter compact manifold, Ai is the four dimensional Minkowski space and both 
A4 and S directions are warped as in the Randall Sundnim model, with two brancs of codimension one sitting at the 
orbifold fixed points. Aside from the usual negative cosmological constant, a bulk sigma model scalar field theory is 
used as the source of gravity in the cases of a curved internal manifold S. We have identified the relevant moduli fields 
characterizing the background, and found the classical action in the moduli approximation (as well as the coupling of 
the moduli to matter fields sitting on the branes). 

We have computed the contribution to the one-loop effective action from generic bulk scalar fields at lowest order 
(i.e. the Casimir energy). The computation, similar to the one for the RS model, is technically more complex since 
there are KK modes propagating along S, resulting in a dependence of the Kaluza-Klein masses on the eigenvalues of 
the Klein-Gordon operator on S. However, for the specific choice of space-time we made, where the warp factors for 
S and for the Minkowski factor A4 are the same, the physical KK masses split as in the usual factorisable geometries. 
Using the Mittag-Leffler expansion for the generalized function we were able to express the Casimir energy in 
terms of heat-kernel coefficients of the internal space S, so that the presence of each divergence is dictated by a 
certain heat-kernel coefficient. This simplifies the renormahzation of the result. An interesting nontrivial check of 
our result is the fact that the RS divergence (which is lower dimensional in this model and which appears as the 
contribution of the S zero mode) cancels out in the final result, as it should, once all contributions are added. We 
renormalized the effective potential by subtracting suitable counter-terms proportional to a number of boundary or 
bulk local operators. Since we work in dimensional regularization, the subtraction is performed in the regularized 
space, with (4 — 2e) + D2 + 1 dimensions. As a result, there is a mismatch in the powers of the moduli appearing 
in the divergent terms, and a number of logarithmic terms (in the moduli) appear in the renormalized expression for 
the effective potential. 

As an application, we proposed a scenario where SUSY is broken at a scale just below the fundamental cutoff M. 
This makes the curvature scale of the background to be several orders of magnitude below M. As a result, a large 
hierarchy is generated by a combination of redshift [6] and a large volume effects [3] . The key point for the latter to be 
efficient (in spite of having codimension one branes) is that the size of the internal manifold S (present in the bulk and 
on the brane) grows as one moves away from the TeV brane, where matter lives. Therefore, this behaves effectively 
as a brane with a small S extra space, attached to which there is a large S space where only gravity propagates. 



The same thing cannot happen for the brane tension terms, since the corresponding coefficient is Co(Pe) = 1 always. 
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As for the stabilization, wc find that, gcncrically, the potential induced by bulk fields can generate sizeable masses 
for the moduli compatible with a large hierarchy with no need of fine tuning if S is flat. If it is curved, the effective 
potential can naturally stabilize the moduli but without a large hierarchy. 

In the model we have considered, the size of the internal space S is everywhere small compared with the size of the 
orbifold. Therefore, there is a range of energy scales where the model is effectively five dimensional (this feature is in 
common with the Hofava-Witten model [4,5]). From the five-dimensional point of view, the model contains a dilaton 
field in the bulk, which causes a power-law warp factor in the Einstein frame, a{y) oc y'^, where y is the; proper distance 
along the extra dimension. The power q is related to the number of additional dimensions through q = {D2 + i)/ D2 
[14], which leads to 1 < q < A. Five dimensional models with power-law warp factors were investigated in [14], where 
it was argued that the counterterms at the orbifold fixed points can naturally stabilize the moduli corresponding 
to the positions of the branes. However, a large hierarchy was not expected unless the power is substantially large, 
q ^ 10. This conclusion is consistent with the results of the present paper, which correspond to relatively small q. In 
this case, the large hierarchy can only be stabilized naturally if the internal space is flat. This case is special because 
the only possible counterterms are renormalizations of the higher dimensional brane tensions. 

The above arguments suggest that a large hierarchy may be obtained by considering more general warped models, 
where a larger power exponent q is obtained after reducing to five dimensions. In such cases, the stabilization of a 
hierarchy without fine tuning is expected even if the internal manifold E is curved and all sorts of comitertcrms are 
present. As a continuation of the work presented here, we are currently considering the case in which the warp factor 
along S is constant [25,29]. We hope to report on this shortly [43]. 
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In the present appendix we prove Eq. (38). It is well known that for a strictly positive definite laplacian P with 
eigenvalues Ap, the associated zeta function C(s|f ) = Ap* admits a Mittag-Leffler expansion of the form 



where Cp{P) are the Seeley-DeWitt coefficients related to P, with p G N/2 and the function f{s\P) is analytic for 
all finite s. This is a very useful relation since it neatly shows the pole structure of the zeta function, in terms of 
geometrical invariants. 

We need to generalize this equation to operators with one zero eigenvalue (go — Consider a positive semidefinite 
differential operator with eigenvalues and assume that there is one zero eigenvalue. As usual in these cases, 
one defines the generalized ( function excluding this eigenvalue (see Eq. (36)), 
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APPENDIX A: MITTAG-LEFFLER EXPANSION 




(Al) 



00 




(A2) 



1=1 



Our main task is to express ({s) in terms of geometrical quantities in the form 




(A3) 
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with p G N/2, for some Cp related to the Seclcy-DcWitt coefficients of Py,- First of all, let us introduce the regulated 
zeta function associated to the operator = + jj? / , 

oo 
1=0 

Now it is trivial to express the function C,{s) in terms of C,^{s), 

as) = lim (C^(5) - , (A5) 

and it is obvious that, understood as this limit, C,{s) is infrared finite for any s, even though Qn^s) is only IR finite for 
Re(s) < 0. By construction, is strictly positive definite, so Cm(*) admits the expansion 



f °° 



^'^^^ (A6) 



where the function g{^, s) is analytic and the Seeley-DeWitt coefficients Cpd^) depend polynomially on /x. 
From Eqns.(A5,A6) it follows that 



1 f °° 



C(») = ^ Ita. O : JT§^ + ») - > . (AT) 



The next step is to isolate the poles from the last term in the previous formula. We do this expanding r(s) as 

(s + 2p)' 



r(s) = r(i,.) + E^-^, (A8) 

p=0 



where r(z, s) is the incomplete gamma function, and the coefficients of the expansion are given by 62^ = (— 
Using (A8) we have: 

00 00 

r(s)M-'« = J2 b2p^^ + E hpif,, s) + r(l, (A9) 

p=0 * p=0 



where we have defined 



Eq. (A9) allows us to write 



— 2s 4p 

hp(p.s) = h / . (AlO) 



where 



/(m, s) = s) - E ^(/x, s) - r(l, s)^x-^', (A12) 

p=0 

and the modified coefficients Cp(^) are then given by 

C'r.,/2+2p(M) = Cd^/2+2p{IJ) - 62p/U^^. (A13) 



'^^Note that this mass is fictitious and has nothing to neither with the physical bulk mass m (16) nor with the renormalization 
scale. 
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Wc note that the only coefRcients which are modified are CB2/21 Cd2/2+i,Cd2/2+2, Taking the Umit ^ 0, we 

obtain the main result of this appendix 



with 



Cp = lim Cj,{ii) = Cp(0) - 5^p,D2 , (A15) 
/ij— >o 



where Sp^pi is the Kronccker delta, the limit of Cp{ji) can be taken because they are polynomials in /i, and the function 
/(s) = /(0,s) is analytic and finite by construction, although it can be explicitly checked from (A12). 

In conclusion, the result (A14,A15) implies that for a laplacian with one zero eigenvalue, there also exists a Mittag- 
Leffer-like expansion for the 'primed' zeta function (A2,36), changing only the Seeley-DeWitt coefficient by 
~ 1- That is why we can consider Eq. (A6) valid in general (with either (70 = 1 or 0), replacing Cyj^ii by 
Cd.212 - 90- 

It is now easy to expand around s = p and a simple calculation gives 

mas) 1,^^ = ^^2i2zE + + o{is pn (aig) 

s p 



with 



APPENDIX B: (D + 1) ^ 5 REDUCTION 



In section VI, wc have argued that there exists a range of energies where the theory is effectively 5 dimensional, as 

illustrated in Fig. 2. 

In this appendix we show the dimensional reduction procedure from {D + 1) dimensions down to 5 dimensions, 
which allows contact with the language of [14]. 

The reduction from the higher dimensional theory (2) to 5 dimensions is performed by the compactification on the 
internal manifold S. This amounts to keeping only the E-zero modes of the fields defined in (_D + 1) dimensions. 

In this section we denote collectively the four dimensional Minkowski coordinates and the orbifold by x°' . 

For the sake of simplicity, we shall consider only the breathing mode of S in the internal components of the metric. 
As well, wc shall freeze the {a,i} components (the graviphotons) to zero. 

Thus, the ansatz for the metric that we shall adopt depends on the internal coordinates AT* only through the 
background geometry on E, and on x" through the five dimensional graviton , and a dilaton cr, 



ds^ = g^{3f)dx"dT^' + R^e^''^'=''^j,jdX'dX^ . 

As for the sigma model scalars, we shall also freeze them to their value in the background, 0" = (/'"(A*). 
The action (2) corresponding to this ansatz is 



(Bl) 



J d^x^e^^- (7^(5) - D2{D2 - l){da)l^) + a} 



(B2) 



where gli^! denote the metrics on the branes induced by and we have performed the X integration. We can 



rewrite this action in the (5 dimensional) Einstein frame, given by 

55 = - Mi J d^x^ S^n, + ^{8^)1 + Ase^^l 



^ dab ^ 



.c4>/2 



(B3) 
(B4) 
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where 

^ 31)2 + 3' 

the canonical scalar field is = — (2£)2/3c)c7, g^v are the metrics on the branes induced by g^^, the 5 dimensional 
Planck mass is given by Af| = v^R'^^M^-'^, A5 = M'^-^h and r5± = v^R^^t±. 

The action (B3) coincides with the 5 dimensional scalar-tensor model considered in [14]. It was found there that 
this model has a solution with a power-law warp factor of the form 

ds^. = a.l{z) [dz^ + i^^t^dx^'dx") , 
Mz) = -V(^M+T)Hz/zo) with a^{z) = {z/zof (B5) 

with p = 2/(3c2 - 2) = -(L>2 + 3)/3. " 

The brane operators induced by quantum effects on this background are given by positive powers of the extrinsic 
curvature scale (see e.g. [14]) IC,^± = P/z±a^± = Pz^^^~^^\ 

/r / \ (4— ")/9 1 r 
d'x^Kl^ = I d'x (^j ^ a I d^X ,4+(4/3)Z..-(V3)i..^ (Bg) 

where n = 1, 2, 3 . . . , we used that the conformal coordinate z = e'^" is the same in 5 and in [D + 1) dimensions, and 
a± = l/kz±. On the other hand, wc have seen in section V that the operators generated by the effective potential 
due to bulk fields in the model (2) are of the form 

I d^x^n^ cx I rf4a;r4+^=-2^ (B7) 

where TZ± is the intrinsic curvature computed with the induced metrics on each brane, g^^. Here = 0, 1, . . . , [{D + 
l)/2], and [ ] denotes the integer part. Now we can identify that these operators correspond to a number of powers 
of the extrinsic curvature operator (B6) given by 

n=^N + l. 

J->2 

We note that all the induced operators can be cast as powers of the extrinsic curvature for D2 = 1,2,3 and 6 only, 
having in the D2 = 6 case a one-to-one correspondence. For any other value of D2, there exist higher dimensional 
local operators that arc not simply powers of /Cg±, but of some power of e"^ in the 5 dimensional effective theory (B3). 

As well, Ref. [14] raised the question that the path integral measure of a bulk scalar field in the effective five 
dimensional theory (B3) quantized on the warped vacuum configuration (B5) is ambiguously defined. The nontrivial 
profile of the scalar (j) permits to define many different conformal frames, all of them equivalent at the classical level. 
However, the path integral measure can be defined covariantly with respect to any of them. It turns out that the 
term proportional to 

In z+ In Z- 



4 



in the potential depends on this choice. Several arguments can be given in favor of possible 'preferred' frames. For 
instance, with a measure covariant with respect to the 5 dimensional Einstein frame metric gf^^, this term is present. 

But if one chooses covariance with respect to g^jj , there is no such term. However, in the model presented here, there 
is no ambiguity in the choice of the measure since in the D + 1 dimensional theory there is no scalar with nontrivial 
profile along the orbifold. In the computation presented here, the choice of the measure shows up (see [14]) when we 
subtract the divergences Eq. (60), covariant precisely with respect to the higher dimensional Einstein frame metric 

S'mjv^'- ^ ^ result, when we take into account both the 5 dimensional modes (the E KK zero mode) together with 
the D + 1 dimensional ones (the KK modes excited along the S as well), we have found that there is a remaining 



'^'*In terms of the proper coordinate (in the 5 dimensional Einstein frame) oc 2;'^+^ , a,^{y,^) = (jJe/VoY with q = 2/3c^ = 

(L»2+3)/£)2. 
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contribution of this form, sec Eqns. (61,55,56). Anyhow, it should be noted that these Coleman- Weinberg-Hke terms 
do not play a very relevant role in stabilizing of the moduli. 
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